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Abstract 

We continue our study of the gravitational collapse of spherically symmetric 
skyrmions. For certain families of initial data, we find the discretely self- 
similar Type II critical transition characterized by the mass scaling exponent 
7 ~ 0.20 and the echoing period A ~ 0.74. We argue that the coincidence 
of these critical exponents with those found previously in the Einstein- Yang- 
Mills model is not accidental but, in fact, the two models belong to the same 

universality class. 
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This is a second paper in the series devoted to the study of critical gravitational collapse 
in the Einstein-Skyrme (ES) model. The main motivation for using this model is an attempt 
to understand the role of different length scales and stationary solutions in the dynamics of 
Einstein's equations at the treshold of black hole formation. In the first paper [[!]] we showed 
that the presence of sphaleron solutions gives rise to Type I critical behavior for certain 
families of initial data. In particular, we found a new kind of first order phase transition in 
which subcritical data relax to the static gravitating skyrmion. Here, we focus our attention 
on Type II critical behavior. After providing the numerical evidence for the existence of 
such transition, we discuss the universality of critical behavior with respect to intrinsic 
scales present in the model. We show that the two length scales become gradually irrelevant 
on small spacetime scales near criticality and consequently the near-critical solutions can be 
asymptotically self-similar. Interestingly enough, the critical exponents characterizing the 
Type II transition, the mass-scaling exponent 7 ~ 0.20 and the echoing period A pa 0.74, 
agree (within the numerical error limits) with those previously found in the Einstein- Yang- 
Mills (EYM) model JJ]. We argue that this fact is a reflection of equivalence of Type II 
critical transitions between the EYM and the ES models, that is the critical solutions in 
both models are asymptotically identical. To our knowledge, this is the first example of 
universality of critical collapse across two physically fundamentally different systems (the 
universality classes observed previously in collapse simulations comprised models differing 
only by potential terms in a lagrangian). 

For this paper to be self-contained let us briefly recall the setup for the spherically 
symmetric ES model we used in Adopting the polar time slicing and the areal radial 
coordinate, the spacetime metric can be written as 

ds 2 = -N(r, t)e- 25{r ^dt 2 + N-\r, t)dr 2 + rW. (1) 

As a matter source we take the SU(2)-valued scalar field U(x) (called the chiral field) and 
assume the hedgehog ansatz U = exp(ia ■ rF(r, £)), where a is the vector of Pauli matrices. 
The components of stress-energy tensor T a \, expressed in the orthonormal frame determined 
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by the metric (1) are 



T a = £ { NF* + N-^ + ^(f> + ^), (2) 



Tn = ^NF° + N-^)-i^(f + °m, (3) 



u 



T 01 = V ^FF\ (4) 

where overdots and primes denote d/dt and d/dr respectively, and u = f 2 r 2 + \ sin 2 F. 
The two coupling constants f 2 and e 2 have dimensions: [f 2 ] = ML^ 1 and [e 2 ] = M _1 L _1 
(we use units in which c = 1). In order to write the evolution equations in the first order 
form, we define an auxilary variable P = ue 5 N~ 1 F. Then, the full set of ES equations is 

F = e" 5 N-, (5) 
u 



P = (e~ & NuF')'-sin(2F)e 



(6) 



1 - N 

N' = SttGVToo, (7) 

r 

N = -8nGre- s NT 01 , (8) 
5' = -4nGrN- 1 (T 00 + T ll ). (9) 

In order to make an identification of certain terms in the equations easier, all the coupling 
constants are displayed, however we remind that, apart from the overall scale, solutions 
depend on these constants only through the dimensionless parameter a = 4irGf 2 . We solve 
Eqs. (5-9) for regular asymptotically flat initial data. The condition of regularity at the 
center N(r,t) = 1 + 0(r 2 ) is ensured by the boundary condition F(r,t) = 0(r) for r — > 0. 
The asymptotic flatness of initial data N(r, 0) = 1 + 0(l/r) for r — > oo is guaranteed by 
the initial condition F(r, 0) = Bit + 0(l/r 2 ), where the integer B, usually referred to as 
the baryon number, is the topological degree of the chiral field. Since the baryon number is 
dynamically preserved, the Cauchy problem falls into infinitely many superselection sectors 
labelled by B - here we concentrate on the B = sector. A typical one-parameter family of 



initial data in this class, interpolating between black-hole and no-black-hole spacetimes, is 
an initially ingoing "Gaussian" 



where one of the parameters A, s, or r is varied (hereafter this parameter is denoted by 
p) while the others are fixed. As usually, the critical value p* is located by performing 
a bisecting search in p. In order to get into close proximity of p* we have implemented 
an adaptive mesh refinement algorithm. This code was essential in probing the critical 
region with sufficient resolution. Our numerical results demonstrate the existence of Type 
II critical transition with its two main characteristic features, first observed by Choptuik in 
the massless scalar field collapse ||, namely: 

• Mass scaling: For supercritical data, the final black hole mass scales as tubh ~ C\p — 
p*\ J with the exponent 7 ~ 0.20. As shown in Fig. 1 this scaling law holds over two 
orders of magnitude of mass. 

• Echoing: For near-critical data, the solutions approach (for sufficiently small r) a 
certain universal intermediate attractor which is discretely self-similar with the echoing 
period A 0.74. This is illustrated in Fig. 2. 

We find that the critical exponents, 7 and A, are universal not only with respect to 
initial data, but also with respect to the parameter a. The universality of Type II critical 
collapse with respect to a dimensionless parameter is a rare phenomenon which can occur 
only if the parameter enters evolution equations through terms which become "irrelevant" 
near criticality (see the discussion of this issue in Gundlach's review ||). Now, we would 
like to show that this is exactly what happens for a in the ES equations. Along the way, 
we will see that discrete self- similarity is compatible with the presence of two length scales 
in the model. Our basic assumption is that F(r, t) / yjf is an echoing quantity (scaling 
variable). This assumption is not only justified empirically (see Fig. 2) but, as follows from 
a simple dimensional analysis of Einstein's equations (7-9), it seems to be the only possibility 




(10) 
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compatible with the discrete self-similarity of metric functions N and 5. We introduce a 
unit of length L = yAwG / e and dimensionless coordinates 

r = h,(*^*) and £ = In (^) , (11) 

where t* is the accumulation time of the infinite number of echos. We also define new 
dimensionless variables 

$(r, = Z (t, = v^L F'(r, f), n(r, = VrLe 5 iV _1 F(r, t). (12) 



By assumption $, Z, and II are the scaling variables, that is they are asymptotically periodic 
in r: $(r, £) ~ $(r + A, £) etc. for large negative r and fixed £. Rewritting Eqs.(5-9) in 
these new variables, we find that a is always multiplied by e T , and the only other terms 
depending explicitly on r appear through the combination 

X = e -^+0 sill ( e !(-+0$). (13) 

For example, the hamiltonian constraint (7) takes the form 

1 _ n - ^ = iV(ae T+ « + 2X 2 )(Z 2 + IT 2 ) + X 2 (2ae r+ « + X 2 ). (14) 

In the limit r — > — oo, the terms containing a become negligible ("irrelevant" in the language 
of renormalization group theory), and therefore the critical behavior does not depend on a. 
In the same limit, X — > $, so the equations become asymptotically autonomous in r, and 
thereby scale invariant. The equivalence of Type II critical behavior between ES and EYM 
models is an immediate consequence of the universality with respect to a, because, as we 
pointed out in ||, for a = Eqs.(5-9) reduce (after the substitution w = cosF) to the EYM 
equations. This means that the EYM critical solution constructed by Gundlach is valid 
(to the leading order) in the ES case as well, and hence the critical exponents in both models 
are the same. As mentioned above, this theoretical retrodiction is confirmed numerically. 

We conclude with two remarks concerning possible extensions of the research presented 
here and in |Q. First, we would like to point out that by setting all terms in Eqs.(5-9) 



containing the coupling constant e 2 (the Skyrme terms) to zero, one obtains the cx-model 
coupled to gravity. This model is scale invariant, so it does not admit a Type I transition 
but a Type II transition is expected to exist. If so, the natural scaling variables will be 

$(r, = F(r, t), Z(t, = rF'(r, t), ft(r, = re 5 iV- 1 F(r, t). (15) 

It is easy to check that in this case the terms proportional to a are not irrelevant. For 
example the analogue of Eq. (14) is 

1-N- — = aN(Z 2 + fl 2 ) + 2a sin 2 ®. (16) 

Therefore, the critical solution, and eo ipso the critical exponents, are anticipated to depend 
strongly on a. We have been informed that the group of researchers led by Peter Aichelburg 
is in the process of investigating this and related problems . If their results confirm our 
expectation, the universality with respect to a in the ES model could be interpreted as 
another nonperturbative effect of the Skyrme correction to the u-model. 

Second, we recall that, in contrast to Type II, the Type I critical transition in the 
ES model is manifestly nonuniversal with respect to a because the critical solution (the 
sphaleron) changes with a, in particular, it exists only for sufficiently small a. Thus, for 
large a only Type II behavior is possible, while for small a the two types of critical behavior 
can coexist. In the latter case, one can anticipate the existence of crossover effects at 
the border of basins of attractions of Type I and Type II critical solutions. We leave 
an investigation of these fascinating effects to other researchers who, as we have recently 
accidentaly found out ||, are also interested in the ES model. We would like to emphasize 
that a full description of an extremely rich phenomenology of the ES model is not a per se 
goal our studies - for us this model is only a testing ground for addressing certain issues of 
the dynamics of Einstein's equations in the presence of intrinsic scales and stationary (stable 
and unstable) solutions. 
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FIGURES 
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FIG. 1. The black hole mass scaling. The logarithm of black hole mass Mbh is plotted versus 
the logarithmic distance from criticality In \p — p*\ for supercritical solutions generated from initial 
data (10) for a = 0.02. The power-law fit is indicated by the solid line with the slope 7 « 0.20. 
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FIG. 2. Numerical evidence for echoing. For a marginally critical solution generated from initial 
data (10) for a = 0.02 we plot the profile of <E> = F/y/r versus p = ln(r) for some arbitrary late 
time ti and superimpose the profile of the first echo at time t2 shifted by p — > p + A p . The time 
ti and the radial echoing period A p are chosen to minimize the mean square difference between 
the two profiles. By repeating this calculation for a sequence of pairs (t i , *2 ) , we estimated the 
temporal echoing period A r from the slope of the line £2 = t*(l — e _Ar ) + e -Ar £i, confirming the 
expectation that A p = A T . 
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